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Abstract. Due to the stochastic nature of biochemical processes, the copy number of 
any given type of molecule inside a living cell often exhibits large temporal fluctuations. 
Here, we develop analytic methods to investigate how the noise arising from a bursting 
input is reshaped by a transport reaction which is either linear or of the Michaelis- 
Menten type. A slow transport rate smoothes out fluctuations at the output end and 
minimizes the impact of bursting on the downstream cellular activities. In the context 
of gene expression in eukaryotic cells, our results indicate that transcriptional bursting 
can be substantially attenuated by the transport of mRNA from nucleus to cytoplasm. 
Saturation of the transport mediators or nuclear pores contributes further to the noise 
reduction. We suggest that the mRNA transport should be taken into account in the 
interpretation of relevant experimental data on transcriptional bursting. 
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1. Introduction 

Molecular binding and chemical modifications underlying intracellular processes are 
intrinsically stochastic. They give rise to temporal fluctuations and cell-to-cell variations 
in the number of molecules of any given type, mask genuine signals and responses, 
and generally contribute to the phenotypic diversity in a population of genetically 
identical individuals [TJ [2]. Various characteristics of such noise have been under intense 
quantitative study in the past few years [HI Hj- One of the focal points of the discussion 
is how the noise propagates along a biological pathway. It has been shown that in 
cases where the dynamics of the upstream molecules is not affected by the downstream 
processes (e.g., in gene transcription and translation), a "noise addition rule" generally 
applies, i.e., each process in the pathway contributes to the overall noise strength 
in a statistically independent fashionj5j E]. Modifications to this rule in molecular 
circuits with feedbacks or "detection" capabilities as in signalling have been examined 
by Tanase-Nicola and his colleagues [7] . 

In metabolic pathways and other transport processes, however, the upstream 
molecule is passed on to the downstream pool in a modified form. The conservation of 
mass and/or number of molecules in the reaction sets new rules on noise propagation. 
In a recent work, Levine and Hwa considered this problem in the steady state of a 
metabolic network [H]. Their results show that fluctuations in the number of intermediate 
metabolites are generally uncorrelated to each other, upstream or downstream along a 
pathway or across branches. This, of course, does not exclude dynamic correlations 
which can be quite nontrivial in driven processes[9]. The full dynamic description 
of stochastic transport through a network is a very challenging task which has not 
been fully solved even for a linear pathway [10]. A case of interest in the present 
context is busty input, where the molecules to be transported are produced in batches 
separated by long silent intervals. Well-known examples of such behavior include 
the transcriptional and translational burstingJH [11] and vesicular transport [12]. In 
eukaryotic gene expression, a large number of mRNAs are produced over a short period 
of time and followed by a long silent period as a result of the chromatin remodeling. In 
prokaryotes, burst of protein copy number can occur as the result of short lifetime of an 
mRNA transcript. Nutrient uptake in endocytosis can also be viewed as a burst event: 
each endocytic vesicle transports and releases a large number of extracellular molecules 
to the target site. 

In this paper, we examine the attenuation of bursting noise in a two-compartment 
model with a stochastic transport channel. The model can be viewed as an improvement 
over the one-compartment model used previously to analyze the mRNA copy number 
fluctuations in single- molecule experiments carried out by Raj et a/.(T3] on mammalian 
cell gene expression. The experiments show that the number of mRNA transcripts 
produced in a single burst event ranges from a few tens to hundreds. It was argued 
that such large bursts, if unattenuated, could harm the progression of normal cellular 
activity[H] due to their large perturbations to the cytoplasmic mRNA and protein levels. 
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A suggestion was made by Raj et al. that the latter could be avoided if the overall 
protein copy number is kept high by a low protein degradation rate. Our analysis here 
points to a second possibility: a slow nuclear mRNA processing and export process 
can also attenuate mRNA bursting and minimize its impact on the downstream protein 
population. 

Supporting this view, an earlier kinetic study of mammalian cell mRNA splicing and 
nuclear transport has shown that the nuclear dwelling (or retention) time of an mRNA 
molecule can be comparable to its lifetime in the cytoplasm [15]. This is evidenced in 
the time required to reach the respective steady state levels for the mRNAs residing in 
the nucleus and in the cytoplasm, which were measured separately in the experiment. 
Consistent with this observation, both studies also concluded that, on average, about 10- 
40% of mRNA are retained in the nucleus. In this respect, the one-compartment model 
of Raj et al, which considers only the total mRNA copy number in a cell, overestimates 
fluctuations in the actual number of mRNAs in the cytoplasm that participate in the 
translation process. 

To establish a reasonable model, let us first examine the typical fate of a single 
mRNA: the mRNA is synthesized in bursts at the transcription site and almost 
simultaneously processed into an mRNA-protein (mRNP) complex [T6l ITT] ; the mRNP 
complex diffuses inside the nucleus [HI [191 120], eventually reaches one of the nuclear 
pores and exits with the help of export mediators |2 11 122] : the mRNA degrades in the 
cytoplasm. Therefore, our model includes three processes: the transcriptional bursting 
in the nucleus, the mRNA transport, and the mRNA decay in the cytoplasm. We 
shall assume that the mRNA transport to the cytoplasm is much slower than the 
diffusion process, so that the spatial inhomogeneity of the mRNA molecules in the 
two compartments can be ignored[23j. 

Both linear and nonlinear transport are considered. In linear transport, the export 
mediators and nuclear pores are abundant as compared to the transported mRNA. 
Consequently, the export events are independent of each other. Nonlinear transport 
corresponds to a situation where queuing of nuclear mRNA takes place due to a limited 
number of transport mediators or nuclear pores. Its mathematical description is identical 
to that of the Michaelis-Menten (MM) reaction in enzyme kinetics. 

The stochastic bursting and transport processes can be described by the chemical 
master equation governing the time-dependent distribution of mRNA copy numbers in 
the two compartments. In the case of linear transport, exact expressions for the copy 
number fluctuations in a steady-state situation are obtained. The MM transport in 
the weak noise limit can be treated using the linear noise approximation (LNA) [24] . A 
novel independent burst approximation (IBA) is introduced to treat the MM transport 
in the strong fluctuations regime. We also perform stochastic simulations to check the 
accuracy of the analytic expressions. 

The main results of our calculation are summarized in Sec. 3. The noise strength 
of cytoplasmic mRNA is expressed in terms of the average burst size and the ratio 
of the mean nuclear and cytoplasmic mRNA copy numbers, which are measurable in 
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experiments. Quite generally, the extent of burst attenuation is governed by the rate of 
transport. The slower the mRNA transport, the smaller is the noise in the cytoplasmic 
mRNA number. In the case of the MM transport, the saturation effect of transport 
mediators or nuclear pores further reduces mRNA copy number fluctuations in the 
cytoplasm. Based on these findings, we suggest a revision of the parameters estimated 
by Raj et al. for the bursting and decay dynamics of mRNA in their experiments. 



2. Methods 



The transport event involves two compartments, usually with different volumes which 
result in different concentrations even when the number of molecules is the same. To 
step aside this problem and to be more clear, we measure the amount of mRNA in copy 
number rather than concentration, and define all the reaction rates mesoscopically(by 
scaling the macroscopic counterparts with volumes). Such a treatment (we will not 
explicitly refer to specific units for variables and parameters in the following calculations) 
also facilitates the application of chemical master equations. Then we can start safely 
by introducing the methods for the simplest case: linear transport. 



2.1. Linear transport 

The burst of mRNA arises from the transitions between active gene A and inactive 
gene /, while each active gene produces an amount of nuclear mRNA M n . To faithfully 
describe the bursting, we write the relevant reactions, as often done in many previous 
works [H [13], as: 

I^A, A-Ul, A^A + M n . (1) 

Here A and 7 are the rates of gene activation and inactivation respectively, and [i is the 
transcription rate. Using M c to denote the cytoplasmic mRNA, we define the linear 
transport and cytoplasmic decay of mRNA as: 

M n M c , M c 0. (2) 

Here k is the transport rate, and 6 is the degradation rate of cytoplasmic mRNA. The 
decay of mRNA in nucleus is ignored |15j . 

At a given time t, the gene can be in either active or inactive state, with the 
probabilities given by P A (m n ,m c ,t) and Pj(m n ,m c ,t), respectively. The probabilities 
depend on the copy number of mRNA in nucleus m n and in cytoplasm m c . The time 
evolution of the two probabilities are governed by the chemical master equations: 
dP A (m Ql m c ,t) 

— = AP/(m n , m c , t) - >yP A {m n , m c , t) 

at 

+ /i^:" 1 - l)P A (m n , m c , t) 

+ kie^' 1 - l)m n P A (m n ,m c ,t) 

+ 5(s c - l)m c P A (m n ,m c ,t), (3) 



Attenuation of transcriptional bursting in mRNA transport 5 
dP I (m n ,m c ,t) 

— = 'jP A (m n ,m c ,t) - AP/(m n ,m c ,t) 

+ k{e n e^ - l)m n P 7 (m n , m c , t) 

+ 5(e c - l)m c P/(m n ,m c ,t), (4) 

where e is the step operator defined by its effect on arbitrary functions of m n and m c : 
^nVOn, m c , t) = f(m n ± 1, m c , t) and e^ 1 f(m n , m c , t) = /(m n , m c ± 1, i). 

The first moments or mean values of m n and m c can be obtained by multiplying 
equations (j3J) and (TjJ by m n and m c in turn, and summing over all m n , m c and gene 
states: 

d(m n ) 

= up A - k(m n ), (5) 



dt 
d(m c ) 



k{m n ) - 5(m c ). (6) 



dt 

Here (■) denotes average over the distribution, and pAif) = J2 mn m c -Pa(^ih^c^) is the 
probability that the gene is in the active state. The following equation is easily seen 
from the gene activation dynamics (CQ): 

dp A /dt = A(l - p A ) - ipa- (7) 

Equations ©-© ar e equivalent to the macroscopic rate equations given by the mass- 
action law due to the linearity of the process. Setting the right-hand-side of these 
equations to zero, we obtain the steady-state relations for the average mRNA flux: 
J = t l P*A = /"V + 7) = k(m n ) = 8(m c ). 

In the following discussion we will focus on the burst limit where fi and 7 are 
significantly larger than all other reaction rates. The steady-state probability p* A 
goes to zero but the mRNA synthesis rate J = fip* A ~ X(^/j) remains finite. Gene 
activation in this case follows a Poisson process at a rate A, but the number of mRNA 
copies b produced in each burst event is a random variable that satisfies the geometric 
distribution: G(b) = (fi/^f) b (l+fi/^)~ b ~ 1 [25]. In terms of the mean mRNA copy number 
produced in each burst, (b) = yu/7, we have 

J = X(b) = k(m n ) = 5(m c ). (8) 

A similar procedure as above yields the second moments of mRNA copy numbers in the 
steady state: 

«) =(m n ) 2 + ((6) + l)(m n ), (9) 

(m n m c ) = (m n )(m c ) + (b) , (10) 

(m n ) + (m c ) 

K 2 ) = K) 2 + ((b) + i)K) - (6) . (ii) 

(m. n ) + (m c ) 

It is customary to measure temporal variations of population size in a stationary 
process using the noise strength (also known as the Fano factor), defined as the variance 
over average [5J [TT] : 
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(m c ) 



(b) + l- (b) 



(m n ) + (m c ) ' 



{m n ) 



(13) 



For the Poissonian fluctuation arising from the simplest case where molecules are 
produced one by one with a constant probability and degraded linearly, the noise 
strength is unity [261 [27] . When the synthesis is burst-like and the degradation is still 
linear, the noise strength becomes (b) + 1, which is much larger than the Poissonian 
fluctuation [28]. This is the case for the nuclear mRNA population [equation (fl2i) ]. 
and for the mRNA without transport, as in the prokaryotic cells. Equation (ITBl . on the 
other hand, shows that although the transport event follows a random process, the noise 
strength of the cytoplasmic mRNA that propagates directly to protein noise, is actually 
reduced. The amount of reduction is controlled by the ratio {m n )/((m n ) + (m c )), which 
increases with decreasing transport rate k. 

2.2. Michaelis-Menten-type transport 

Accumulation of the mRNAs in the nucleus may lead to a saturation effect that changes 
the transport dynamics when the number of export mediators or nuclear pores becomes 
limiting. This prompts us to study a more general mechanism of transport that takes 
the transport capacity into account. The resulting process can be cast in the form of 
the well-known Michaelis-Menten model. 

To simplify the discussion, we consider here only one source of constraint, say the 
limited number of one kind of export mediator denoted by E, and treat the rest of the 
export mediators (including the nuclear pores) in the process to be non-rate-limiting. 
Thus, the transport process of mRNA can be described as: 



where k\ and k% are the binding and unbinding rates respectively, and k% is the export 
rate. This is similar to the MM model for enzymatic reaction. 

The analysis presented below is based on the "fast equilibration" approximation, 
in which EM n is treated as a transition state rather than an accumulation point in 
the mRNA export. For this to be true, the lifetime of the complex EM n should be 
significantly shorter than the total nuclear dwelling time, i.e., either the Michaelis 
constant K m = (/C2+&3) / k\ is much greater than one, or k% is much smaller than the other 
two rates. Under this assumption, the EM n population remains in quasiequilibrium with 
the nuclear population m n which varies on a much slower time scale as compared to the 
decomposition time of the complex EM n . 

Due to the small copy number of the nuclear mRNA and the transporter E, we 
distinguish the free M n from the bound ones, and define the total number of nuclear 
mRNA as m n = m n f + c, where m n f is the number of free M Q and c is the number of the 
complexes. The total number of E, including both free and bound ones, is denoted by 
e t . The usual rate equation approach for the complex yields /cim nf (e t — c) — (k 2 + k 3 )c = 



M n + E^±EM 1 



11 



E + M< 



(14) 
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when quasiequilibrium is established. Hence the mean complex number is given by, 

c=e«-^. (15) 
K m + m nf 

On the other hand, through an exact analysis, Levine and Hwa[8] obtained a modified 
expression in the limit k$ — > 0, 

c = e t — , 16 

K + m n 

where K = K m + e t . The two expressions converge to the exact result in both the linear 
regime m n <C K and the saturated regime m n 3> K. They differ only in the crossover 
regime m n ~ K, where no exact result is available in the general case, though either of 
the two can be used as approximate expressions. 

Using equation ffl6|) . we write the transport flux as, 

v(m n ) = k 3 c = — , (17) 

K + m n 

where v max = k 3 e t . Following the fast equilibration assumption, we may now describe 
the mRNA export under the MM kinetics using the reduced description ([2]) with an 
effective transport coefficient k = v m3jX /(K + m n ) that decreases with increasing m n . 
Again, taking the burst limit for the mRNA production, we arrive at the following 
master equation for the joint distribution of m n and m c : 

dP(m^m c ,t) = g XG{b)p{mn _ ^ t) _J- \G(b)P(m nj m c , t) 

b=0 b=0 

+ (ens" 1 - l)v(m Q )P(m n ,m c ,t) 

+ 5{e c - l)m c P(m n ,m c ,t). (18) 

Since the nonlinear function v(m n ) in the MM transport does not allow for the closure 
of equations for the second moments of the distribution, approximate treatment of the 
master equation is necessary. 

The two limiting situations (m n ) ^> (b) and (m n ) <^ (b) call for different 
considerations. In the former case, the contribution of each burst event on the total 
nuclear mRNA population is small, so that perturbative treatment around the average 
(m n ) is appropriate. The latter case, however, corresponds to the situation where the 
nuclear mRNA from each burst event is essentially cleared before the next one arrives. 
The two cases are treated separately below. 

2.2.1. Linear approximation We first consider the weak fluctuation case (m n ) ^> (b). 
A general scheme to perform the noise calculations is the van Kampen's f2-expansion 
whose lowest order terms reproduce the macroscopic rate equations and the next order 
terms yield a linear Fokker-Planck equation (FPE), which is often called the linear noise 
approximation [21]. In Appendix A we derive the noise strengths under the LNA. Here, 
we outline a more direct yet equivalent way to obtain the results. 
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The approximation we introduce is to replace (1171) with its linear expansion at 
m D = (m n ): 

v ~ k eS (rn n + m ). (19) 

Here k c s = v milx K/(K + {m n )) 2 and m = (m n ) 2 /K. Substituting (TT9T) into (|T8|) . we 
may compute moments of the distribution P(m n ,m c ,t) in the same way as in Sec. 
2.1. In fact, equations (1H1)- (1TT1) remain valid if we make the substitution k —>■ k e s and 
m n — > m Q + m . After rearranging the terms, we obtain: 



(m n ) K 



!)«&> + 1), (20) 



w 

)_ 

K+(m n ) 



= (6) + 1 - (6) y . n/ . (21) 



2.2.2. Independent burst approximation In the case (m n ) (6), the mRNA produced 
in a given burst has sufficient time to exit the nucleus before the next burst arrives. It 
is then appropriate to consider the independent burst approximation, where individual 
burst events contribute additively to m n (t) and m c (t): 

= £)£(M-*i), (22) 
m c (t) = J^viKt-U). (23) 

Here 6, is the size of the ith burst which takes place at ij, and £(b, t) and 77(6, t) are the 
number of mRNAs in the nucleus and in the cytoplasm, respectively, generated by a 
single burst of size b at t — 0. 

Three independent stochastic processes contribute to the statistical properties of 
the time series m n (t) and m c (t): i) the time of the burst events tj, which we assume 
to be Poisson at a rate A; ii) the size 6; of individual bursts which follow the geometric 
distribution G(b); and iii) the stochastic nature of the MM transport and mRNA decay 
in the small copy number regime. In the following discussion we shall focus on the 
noise effects due to processes i) and ii), while neglecting stochasticity in iii). The 
latter approximation is justified by noting that the most significant contributions to 
the quantities computed below are from the period when £(&, t) and rj(b, t) are large and 
their relative fluctuations are small. Denoting by x(b, t) = (£(b,t)) and y(b,t) = (rj(b,t)), 
we obtain the following expressions for the moments: 

POD 

(m n ) = J2 G ( b ) / x(b,t)Xdt, (24) 







m l) = ( m n) 2 + 



poo 

VG(t) / x 2 (b,t)Xdt, (25) 
, Jo 



b 

f 00 



POO 

(m c ) = VG(t) / y(b,t)Xdt, (26) 
Jo 
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b J ° 

POO 

(m n m c ) = (m n )(m c ) + > G{b) / x(b,t)y(b,t)Xdt. (28) 

Under the MM transport ( JT71) . the dynamical equations for x(b,t) and y(b,t) are given 
by: 

~77 — ~~ ^maxTT^ , (29) 

at A + x 

| = " ^ 0°) 

with the initial condition 0) = b and ?/(&, 0) = 0. 

As shown in Appendix B, the sum and integrals in equations (124j) and (1251) can be 
worked out exactly to give: 

(m n ) = -^(b}((b} + K+ 1 -), (31) 

"max ' * ' 

<4 n =(m2)-(m n )2 = ^( 6 )[2(6) 2 + (A + 2)(6) + | + il. (32) 



2 3 



Hence, 



(m n ) " {b} + 2 + A + \ + (b) ■ (33) 

The mean value of m c can be obtained from flux balance, i.e., A(6) = 5(m c ) or 
(m c ) = (\/5)(b). The calculation of (m 2 c ) is a bit more involved which we relegate to 
Appendix B. Assuming 5 <C v max , i.e., the decay time of an mRNA molecule is much 
longer than the fastest release time of one mRNA to the cytoplasm, we may write the 
result in the form: 

/ (b) + l)^(u,w), (34) 



K) V w 2, 

where u = (b)5/v max and w = KS/v max . The function \1/ is given by 



1 pi gU)ln(x/xi) 



ty(u,w) = 2 dx dxij- -. tt-. (35) 

Jo L [1 + M(X!-X)] 3 

Since the integrand is less than 1, we have *&(u, w) < 1. 

We have not managed to find a closed form expression for \&(it, w), but the integral 
can be worked out in the two limiting cases: i) u — 0,^(0, w) = 1/(1 + w); and ii) 
w = 0,^(u,0) = 1/(1 + u). These two expressions also set upper bounds for ty(u,w) 
in general. An approximate expression that is consistent with the two limits and also 
verified by numerical integration of fl35l) is given by: 

*(u,w)~ 1 ■ (36) 
1 + u + w 
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In terms of (m n ) and (m c ) and with the help of (131 p and the flux balance condition, 
equation flMl) can be rewritten as: 



a 



2 



(m c ) l + K +<b) 

X ' ^ (m c ) K+(b) + \ 



(37) 



3. Results and discussion 

Let us first summarize the analytical results derived in Section 2 when the average burst 
size (b) 3> 1. In general, noise strength of mRNA copy number in the two compartments 
can be expressed in the form, 

9 

^r = «(6) + l, (38) 

= n } h \ ll r + 1. (39) 

(m c ) 1 +fj(m n )/(m c ) 
For the linear model, we have en = (3 — 1. In the MM case where the transporter may 
become the bottleneck in the process, a = (3 = 1 + (m n ) / K if the slow transport leads to 
the nuclear accumulation of mRNA, i.e., (m n ) 3> (b). In the opposite limit (m n ) <C (b), 
where there is nuclear clearance between successive burst events, a ~ 1 + (b) / (K + (&)) 
and /3 ~ 1. 

The general trend of noise attenuation on the cytoplasmic mRNA due to delay in 
nuclear transport is now clear. Retention of the mRNA inside the nucleus decreases the 
effective burst size and hence the noise strength of the cytoplasmic mRNA. Even at a 
fixed ratio of (m n ) to (m c ), further reduction of the noise is possible in the nonlinear 
MM transport when (m n ) is greater than both the dissociation constant K and burst 
size (b). On the other hand, the independent burst approximation yields a (5 value 
close to one, extending the validity of the linear model when the noise strength of m c is 
considered as a function of the ratio of mean copy numbers (m Q ) / (m c ). 

Fluctuations in the nuclear mRNA level, on the other hand, exhibits a somewhat 
different behavior. The parameter a that characterizes the noise strength of m n reaches 
the minimum value 1 in the linear model, but increases when the saturation effect in 
the MM transport kicks in. Thus queuing results in an enhanced fluctuation upstream 
of the transport channel. 

To check the accuracy of the analytic results under parameter values that broadly 
correspond to the mammalian cell gene expression experiments mentioned above, we 
have carried out simulations of the stochastic MM transport defined by (THl) . following 
the Gillespie's exact algorithm |29j. The unit of time is chosen such that the mRNA 
decay rate 5 = 1. The number of transport channels is set to e t = 10. We fix the 
average mRNA production rate A (b) relative to the mRNA decay rate S such that the 
mean cytoplasmic mRNA copy number (m c ) = X(b)/5 = 40. For easy comparison 
with the experimental measurements and analytic results, we also fix the mean nuclear 
mRNA copy number (m Q ) = 20. With these parameter values, the noise strengths in 
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(a) <b>=4 (b) (b>=20 (c) (b>=80 




Figure 1. Noise strengths (normalized by their corresponding values in the linear 
model) of the nuclear (upper panel) and the cytoplasmic (lower panel) mRNA copy 
numbers against the Michaelis constant K m which controls saturation of transport 
channels. Here e t = 10, (m n ) = 20 and (m c ) = 40. Results are shown for stochastic 
simulation (open circles), LNA (solid line) and IB A (dashed line). The overall noise 
strength is set by the mean burst size (b) which is chosen to be (in units of mRNA 
copy number): 4 [(a) and (d)]; 20 [(b) and (e)]; 80 [(c) and (f)]. See text for the choice 
of other model parameters. 

the linear model are given by o^J (m n ) = (b) + 1 and c^ c / (m c ) = + 1, respectively. 
Simulations are then performed to examine the effect of channel saturation on the noise 
strengths by varying the MM parameters k\ and k?, in such a way that the mean nuclear 
mRNA copy number stays at the value set above. The unbinding rate of the mRNA- 
transporter complex EM n is fixed at a low value k<i — 0.1. 

Figure [1] shows a comparison of simulation data (open circles) and analytic results 
under the LNA (solid line) and the IBA (dashed line), respectively. Plots on the upper 
panel give the noise strength of nuclear mRNA as a function of the Michaelis constant 
Km — (k 2 + k 3 )/ki. Fluctuations in m n grow as the saturation effect becomes more 
prominent on the low K m side. An opposite trend is seen in the fluctuations of the 
cytoplasmic mRNA copy number m c shown in plots on the lower panel. As expected, 
the LNA results (solid line) agree well with the simulation data (circles) in the weak 
burst regime [(a)and (d)], in which case the overall noise strength (as compared to the 
mean mRNA copy number) is weak. On the other hand, the IBA results (dashed line) 
represent a better approximation in the strong burst regime [(c) and (f)][30j. Therefore 
each of the two approximations perform reasonably well in their respective regimes of 
validity, and are complementary to each other. 

Figure [2] shows the actual distribution [(a) and (c)] and sample time course [(b) and 
(d)] of the nuclear and cytoplasmic mRNA copy number, respectively, generated from 
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(a) 



x 10 



(b) 




10 15 20 
Time 



Figure 2. Distributions of mRNA copy numbers in the nucleus (a) and in the 
cytoplasm (c), obtained from 500 000 repeated runs of the Gillespie exact simulation. 
Sample time courses of m n and m c in a single simulation are shown in (b) and (d). 
Parameter values are: A = 2, (b) = 20, k% = 1.56, h-x = 0.1, = 10, e t = 10 and 
S = 1. 



simulations at K m = 6.5. Other model parameters are the same as that of figured] (b) 
and (e), which represents a borderline case for the two approximate treatments. It is 
seen from figure [2] (b) that m n falls below K = 16.5 most of the time, so approximating 
the MM transport by a linear expansion at m n = (m n ) is too crude. On the other hand, 
there are occasional overlaps of the mRNA produced in successive bursts. This has the 
effect of slowing down the mRNA transport than what is assumed in the IBA, leading 
to a somewhat lower m c noise as seen in the left part of figured] (e). 

We have also examined the validity of equation (TiTl) for the mean transport flux. At 
a given number c of complexes EM n , the expected mRNA export flux is k^c. Therefore 
v(m n ) can be obtained in the simulations from the conditional average of c at a given 
m n . As shown in figure [3] equation (fT7|) (dashed line) fits the simulation data (dots) 
very well when K m > (m n ) [(c)], but quite poorly in the regime K m (m n ) [(a)]. 
Surprisingly, the classic MM equation (TT5|) . suitably modified to be considered as a 
function of m n , agrees with the numerical result extremely well for both large and 
small values of K m . Note that equation ffTTl) was derived under the quasi-equilibrium 
approximation k\,k2 which does not hold in the present case. This inaccuracy 

may also contribute to the discrepancy between our analytic results and simulation data 
at small K m in figure [1] 
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Figure 3. Transport flux as calculated from the stochastic simulations (dots), equation 
(fT5|) (solid line) and equation (JT7J) (dashed line). The transport parameters are: (a) 
k 3 = 8.4, K m = 0.1; (b) fc 3 = 14, K m = 24.5; (c) k 3 = 30, K m = 102.5. Other 
parameter values are given by A = 2, (b) = 20, k-2 = 0.1, et = 10 and 5=1. Note that 
(m n ) = 20 and (m c ) = 40 are the same as before. 



4. Conclusions and outlook 

The main conclusion of our work is that the nuclear envelope, which sets a natural 
barrier for the exodus of mature mRNA to the cytoplasm in an eukaryotic cell, can 
significantly attenuate the effect of transcriptional bursting on the downstream protein 
population. The extent of the noise reduction on the cytoplasmic mRNA copy number 
is controlled by the transport efficiency. A high transport rate has a weak effect in noise 
reduction and essentially brings one back to the one-compartment model considered 
previously. On the other hand, a low transport rate turns the nucleus into a buffer 
for the bursting noise, thereby reducing the temporal variation of the mRNA copy 
number in the cytoplasm. This effect is more dramatic in the saturated regime under the 
Michaelis-Menten dynamics where, due to the limited availability of transport channels, 
the mRNA export becomes a Poisson process unaffected by the bursty input. 

Our results can be used to re-estimate the mean mRNA burst size in the experiment 
by Raj and his colleagues. Take the linear transport as an exmaple, the noise strength 
of the total mRNA copy number in the cell can be obtained from equations ([9]), (fTOl) 
and (TIT]): 

- t al T c \ = <*> + 1 + ® u ( Tl ( r\^ - ( 4 °) 

(m n + m c ) {{m n ) + {m c )Y 

Thus the average burst size (b) can be obtained from the measured total mRNA 
fluctuations and the ratio (m n ) / (m c ) of nuclear to cytoplasmic mRNA. If 30% of the 
mRNA accumulate in the nucleus, the burst size should be 83% of that estimated from 
the model without transport. With the help of the two-compartment models, it would 
be interesting to revisit the single-molecule experiment data of Raj et al. to gain a 
more complete view of the role of transport on the characteristics of the mRNA noise 
generated by transcriptional bursting. 

Previous studies suggest that for eukaryotes, such as S. cerevisiae\2\, 
Dictyostelium^yQ and mammalian cells [13], transcriptional bursting is a dominant source 
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of the noise from the internal molecular circuits known as the "internal noise". It is 
interesting to note that the stress-related genes, which demand a fast response and 
an accelerated export rate, are noisier than the essential house-keeping genes such as 
proteasome genes [32]. In our limited transport capacity model described by the MM 
dynamics, the downstream noise is indeed stronger in the linear regime and weaker in 
the saturated regime. The linear regime also allows for a faster response to external 
stimulus as there is no queuing effect. In this respect, one can not help but wonder 
whether special processing and export channels are in place in the nucleus for the fast 
release of stress response genes without congestion. 

Finally, we would like to mention that noise propagation through capacity-limited 
channels with a bursty input is a general phenomenon not limited to the transcriptional 
bursting, and hence the discussions initiated here can be of broader significance. A 
somewhat peculiar behavior associated with the limited-capacity transport is that, 
as the extent of channel saturation increases while maintaining the same average 
transport current, the noise level of the upstream population increases while that of 
the downstream population decreases (see figured]). This counterintuitive phenomenon 
is neither stochastic resonance (SR) nor stochastic focusing (SF) proposed by Paulsson 
et al [26j EZ] • SR is usually related to periodic signal detection, where the signal noise is 
typically external, while SF exploits signal noise to make a gradual response mechanism 
work more like a threshold mechanism. Here, the interesting phenomenon we observe is 
due to a quite different mechanism: the degradation of upstream species and production 
of downstream species share a common reaction (here transport), whose effective order 
at steady state can be tuned from zero to one by a combination of reaction parameters. 
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Appendix A. The fi-expansion 

Equations (|20|) and ([2~T]) can be equivalently derived using the more formal f2-expansion 
which applies when fluctuations of m n and m c are weak. (Q here stands for the system 
volume.) To set the notation straight, the rate equation for the copy number X{ of 
molecule i (i.e., m n or m c in the present case) is given by, 



where Si q is the stoichiometric coefficient of molecule i in reaction q, and V q (which 
depends on the copy numbers in general) is the propensity of reaction q. Stochasticity 
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in the problem is assumed to be due to fluctuations in the propensity V q with a strength 

1/2 

set at V q due to the underlying Poisson process. 

In the large volume limit and for the steady state, the joint distribution of the x»'s, 
which satisfies the Fokker-Planck equation, can be approximated by a gaussian centered 
at the solution to the equation dx./dt = S • V = 0. The width of the distribution is 
parametrized by a covariance matrix C which satisfies the equation [33] : 

AC + CA T + B = 0. (A.2) 



Ajj = (A.3) 



Here 

d(S ■ V) t 

<? 

all evaluated at the steady state. Specializing on our two-component problem, there 
are three reactions: the bursting reaction leading to the production of m n , the MM 
transport reaction for nuclear export, and the mRNA decay in the cytoplasm. Simple 
calculations yield: 



' (K+(m„» 2 U » 



Application of ( 1A.4I) yields: 

B n = £xtH + E ^G{b)b 2 
K + (m n ) ^ 

A + (m n ) 

= 2A((6) + (6) 2 ). (A.6) 

In the last step, we have used the flux-balance condition and (b 2 ) = 2(6) 2 + (b). Other 
matrix elements of B can be obtained as: 

B n = B 21 = -X(b), B 22 = 2X(b). (A.7) 

Substituting the values of A and B into equation ( 1A.2I) gives: 

ok =C n = ^^ +K \ { b) + l), (A.8) 

(m n m c ) - (m n )(m c ) = C 12 = C 21 = (b) K ^ ){m f - , (A.9) 

< =Q a = K)«6) + i)-<6) r . (A.io) 



(A.5) 



K(m c ) 
K+(m n ) 

from which equations fl20l) and (|21J) follow. 
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Appendix B. Integrals in the independent burst approximation 

A convenient way to carry out the integrals in (124l) - (127l) is to convert them into 
integration over x, which decreases monotonically from its initial value b to in a 
single burst event, with the help of (129]) and ( l30l) . Following this procedure, we may 
write, 

dtx(b,t) = - [ b ^r= f dx^^ = —{Kb+-b 2 ). (B.l) 



/q dx / dt Jo f max fmax 2 

, 2/I , f b x 2 dx f b , K + x 
dtx {b,tj — — I - — — - = / xdx 

Jo ax/at Jo f max 

^/rt 2 + -6 3 ). (B.2) 



max 



To perform the averaging over b, we make use of the following results for the geometric 
distribution, 

<6 2 > = <5>(l + 2<6», 

(b 3 ) = (b) (1 + 6(b) + 6(b) 2 ). 

With the help of these results, fl3Tj) are fl32|) are readily obtained. 
The dependence of y on x follows the equation, 

^ = ^Z| i = _i + _*_(! + £),, (B.3) 

which can be integrated to give, 

r b 

y(x) = / dx^-'^ 1 ^ +™ ln (*M), (B.4) 

J X 

where w = SK/v m£LX . 

As a check, let us first consider 

f°° f b K + x 

/ dty = / dx y. (B.5) 

70 JO "max 1 

Using equation (IB .3|) and noting that y(x = b) = y(x = 0) = 0, we obtain, 

f dx^^y = cT 1 f dx(l + ^-)= b/5. (B.6) 

JO ^maxl Jo dX 

Hence, 

K> = \(b)/5, (B.7) 

which is nothing but the conservation law. 
We now consider 

r b dv r b 

dty 2 = / y5~ 1 (-^ + l)dx = 5- 1 ydx. (B.8) 



' dx 

Using equation (1B.4|) and perform the substitution x — > bx, we obtain 



/■OO 1,2 /-l /•! 

/ rfty 2 = _ / dx d Xie &Wtw)(*-*i)+w]B(*/*i). (B.9) 
Jo o Jo A 
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The averaging over b can now be readily carried out. Using the result YlbLo ^ 2 ° ft = 
a(l + a)/(l — a) 3 , we obtain, 

f°° r 1 r 1 p w\n(x/xi) n (-t _i_ n \ 

X>(6) jf *tfM = r*J t * I ^ 1 + {b) ( |_ + o) i . (B.io) 

Here a = (&)e< x_xl W w /(l + (6)). 

To avoid run-away accumulation of mRNAs in the nucleus, we require X(b) < v m3jX . 
Therefore (m c ) < v max /5. Note that (m c ) > 1 automatically implies S/v max to be a 
small quantity In this case, we can approximate a ~ (b)[l + {x — Xi)8 /v max \/{l + (b)). 
Consequently, 

J2G(b) f°<% 2 (M) = mb l + ^ *(u : w), (B.ll) 

where u = {b)5/v majX and ^/(u,w) is given by ( |35l ). 

Finally, the integral in equation ( 1281) can be rewritten in the form, 

poo rb j r2 pb 

/ dtx(b,t)y(b,t) = / a;r 1 (-^ + l)dx= — -r 1 / yrfx. (B.12) 
Jo Jo dx 26 J 

Comparing with ( IB. 81) and using (IB. lip , we obtain, 

VG?(6) /°° dto(6, t)y(6, t) = (6) + ^ [1 - ^(«, (B.13) 
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